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ABSTRACT 

We address the origin of the robust bi-modahty observed in galaxy properties about a 
characteristic stellar mass ^3x lO^^M©. Less massive galaxies tend to be ungrouped 
blue star-forming discs, while more massive galaxies are typically grouped red old-star 
spheroids. Color-magnitude data show a gap between the red and blue sequences, ex- 
tremely red luminous galaxies already at 2:~ 1, a truncation of today's blue sequence 
above i*, and massive starbursts at z ~ 2 — 4. We propose that these features are 
driven by the thermal properties of the inflowing gas and their interplay with the 
clustering and feedback processes, all functions of the dark-matter halo mass and as- 
sociated with a similar characteristic scale. In haloes below a critical shock-heating 
mass Mshock~10^^-MQ, discs are built by cold streams, not heated by a virial shock, 
yielding efficient early star formation. It is regulated by supernova feedback into a 
long sequence of bursts in blue galaxies constrained to a "fundamental line" . Cold 
streams penetrating through hot media in M ^ Mjhock haloes preferentially at z ^ 2 
lead to massive starbursts in L >L, galaxies. At z<2, in M >Mshock haloes hosting 
groups, the gas is heated by a virial shock, and being dilute it becomes vulnerable 
to feedback from energetic sources such as AGNs. This shuts off gas supply and pre- 
vents further star formation, leading by passive evolution to "red-and-dead" massive 
spheroids starting at z^l. A minimum in feedback efhcicncy near Mghock explains the 
observed minimum in AI/L and the qualitative features of the star-formation history. 
The cold flows provide a hint for solving the angular-momentum problem. When these 
processes are incorporated in simulations they recover the main bi-modality features 
and solve other open puzzles. 

Key words: cooling flows — dark matter — galaxies: evolution — galaxies: formation 
— galaxies: haloes — shock waves 



1 INTRODUCTION 

Observations reveal a robust bi-modality in the galaxy 
population, being divided into two classes, the "blue" 
and "red" sequences, at a characteristic stellar mass 
Ms^crit — 3 X W^^Mq. This corresponds to a dark-halo mass 
Merit ~ 10^^ Mq and a virial velocity Krit — 120kms~^ to- 
day. Less massive galaxies tend to be blue, star-forming discs 
residing in the "field" . Their properties are correlated along 
a "fundamental line" of decreasing surface brightness, in- 
ternal velocity and metallicity with decreasing luminosity. 
Galaxies above Merit are dominated by spheroids of red, old 
stars, with high surface brightness and metallicity indepen- 
dent of luminosity. They tend to reside in the high-density 
environments of groups and clusters and they preferentially 
host Active Galactic Nuclei (AGN). 

Current models of galaxy formation have difficulties in 
reproducing this bi-modality and the broad color distribu- 
tion observed. In particular, the extremely red bright el- 



lipticals which start showing up already at 2 ~ 1 are not 
predicted. They require efficient star formation at earlier 
epochs, followed by an effective shut-down of star forma- 
tion in massive galaxies. The observations also reveal very 
blue galaxies in excess of the predictions, indicating repeat- 
ing starbursts over the lifetimes of galaxies. Today's blue 
sequence is non-trivially truncated at the bright end, while 
at 2: ^ 2 there are indications for very luminous starbursts 
in big objects, both posing severe theoretical challenges. 

1.1 The observed bi-modality 

The bi-modality or transition in galaxy properties is ob- 
served in many different ways. We list here the main relevant 
observed features^ which we address in this paper 

^ Quoting only sample references, making no attempt to be com- 
plete. 
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(a) Luminosity functions. Blue galaxies dominate the 
stellar mass f unction below Mg.crit while red gala xies take 
over above it jBaldrv et al.l l2004: Bol l et aljE)03bl in SDSS 
and 2MASS). The transition occurs slightly below L,, the 
characteristic luminosity of the brightest disc galaxies, be- 
yond which the luminosity function drops. 

(b) Color-magnitude. A color bi-modality shows robustly 
in color-magnitude diagrams, where the galaxies are divided 
into a blue sequen ce and a re d sequence separated by a gap. 
In SDSS (Blanto n"eraII 1200 3: Baldrv et al. 2004), the gap 
is at u — r ~ 2. The color distribution is non-trivially broad, 
with the red tip stretching beyond u — r — 2.5 and the blue 
tail reaching well below u — r = 1.0. 

(c) Star-formation rate. The current star-formation rate 
(SFR), and the typical age of the stellar population, show a 
robust bi-modality about Ms,crit. The less massive galaxies 
are dominated by young populations, wh ile the more mas- 
sive galaxies are dominated by old stars (iKauffmann et aD 
l2003l:lMadgwick et all2003l. in SDSS and 2dF), in agreement 
with the color bi-modality. A similar bi-modality is seen in 
the gas-to-stellar mass fraction, which is high in the blue se- 
quence and low in the red sequence, steeply increasing with 
stellar mass below Ms,crit, and only moderately so above it 
jKannaDpan"2004. in SDSS+2MASS). 

(d) Color-magnitude a t z ~ 1. The col or bi-modality is 
similar back to z ~ 1.5 jBell et aljliooi . in C0MB017). 
Extremely red massive galaxies (EROs) exist at the 
bright tip of the red sequence already at z ~ 1 (e.g. 
iMoustakas et al1l2004ll . Very blue small galaxies indicating 
starbursts show in the blue sequence (e .g. lFerguson fc Babull 
ll998l:lFioc fc Rocca-VolmerangeHigogll . 

(e) Massive starbursts at high z. Very luminous and 
massive dusty objects are detected at z ~ 2 — 4, indicating 
an ex cessive activity of star formation in surprisingly big ob- 
jects (IShaplev et all2o'oi : ISmail et alj2002l : IChapman et alJ 
120031.12004 LBG and SCUBA sources). 

(f) Star-formation history. The cosmological history of 
star- formation rate has a broad maximum near z ~ 1 — 2, fol- 
lowed by a sharp drop from 2: ~ 1 to z = fe.g.|M adau e1^dj 
19961 : [Dickinson et a,ijl2003tlHartwickl2004l:lG4avalisco et'all 



20041 : iHeavens et a0l2004^ . StiU, about half the stars in to- 
day's spirals seem to have formed after 2 ~ 1, e.g. in lumi- 
nous infrared galaxies (LIRGs) near Mg^crit iHammer et alJ 
I2OO4I) . Massive galaxies tend t o form their stars ear lier than 
smaller galaxies ( "downsizing" iThomas et alJl2005^ . 

(g) Bulge-to-disc ratio. The galaxy bulge-to-disc ratio 
shows a transition from disc dominance in the blue sequence 
below Ms.crit to spheroid dominance in th e red sequence 
jKauffmann et al.ll2003l : iBlanton et al.ll2003l . in SDSS). 

(h) Environment dependence. The distributions in color 
and SFR depend strongly on the galaxy density in the 
~ IMpc vicinity: the blue and red sequence galaxies tend 
to populate l ow and high density environme n ts respectively 
jHogg et al.l l2003l: IKauffmann et"ai] l2004l: Isianton et alJ 



1200a: iBalogh et alj 120041: iBlanton et al.l |2004 in SDSS). 
The color-environment cor relation i s stro nger than the 
morphology-density relation [Pressleil il980t) . 
(i) Halo mass and HOD. The environment density is cor- 
related with the mass of the host dark-matter (DM) halo, 
where haloes less massive than ~ 10^^ Mq typically host 
one dominant galaxy each while more massive haloes tend 
to host groups and clusters of luminous galaxies, as quanti- 



fied by the Halo Occupation Distribution fHOD.'Ya n et all 
[2003; Abazaiian ct al. 2004; Kravtsov ct al. 2004, in 2dF, 
SDSS and in simulations). The environment dependence 
thus implies that the galaxy properties are correlated with 
the host-halo mass, w ith the bi-modality at Merit 10^^ Mq 
iBlanton e t al."2004\ 

(j) Hot halo gas. Ellipticals of Lb > 10^"'^ show 
a significant exce ss of X-ray flux plausibly associated 
with hot halo gas iCiotti et alJIiggH : iMathews fc Brighentil 



120031) . Inter-galactic X-ray radiation is detected predomi- 
nantly in groups where the brightest galaxy is an elliptical. 
Group properties have a transition near gy ~ 140 km s~ ^ 
jHelsdon fc Ponmanll20o3 : lOsmond fc Ponmanll20o3) . 
(k) Luminosity /mass functions. The stellar-mass 
function has a "knee" near A/s,crit, where the shal- 
low dn/dMs oc M^^ on the faint side turns into an 
exponential drop. In contrast, the dark-halo mass func- 
tion is predicted by the standard ACDM model to be 
dn/dM (X M"^ * everywhere below ~ 10^^ Mq. A match 
at Ms,crit requires a baryonic fraction Ms/M ~ 0.05, 
indicating gas loss, and associating Ms,crit with 
Merit ~ 6 X 10" M0. The halo mass-to-hght function 
has a minimum near Aferit, varying as M/L oc M~'^^^ 
and oc M"*'^''^ below and above it respectively, thus im- 
plying increased suppression of star f ormation away from 
the critical mass on both sides iMari noni fc Hu dsonI 
20021: iBeU. Baugh. Cole. Frenk fc LacevI l2003| : 

" Bell. Mcintosh. Katz fc Weinberd' l2003al : 



Yang. Mo fc van den BoschI I200I in a B-mag sample 



and in SDSS, 2MASS, 2dF). 

(1) Fundamental line vs. plane. A transition is detected 
in the galaxy structural scaling relations near Ms,crit, e.g., 
the surface brightness changes from /is oc M°'^ at lower 
masses to /is ~ const, at the bright end iKauffmann et all 
I2OO3I in SDSS). The correlation below Ms,crit is part of the 
"fundamental line" relating stell ar mass, radius and veloc- 
ity over five decades in Ms (e.g. iDekel fc W oo"2003>l. The 
mean metallicity shows a transition near a similar mass 
scale from Z M? "^ to Z ~ const. llTremonti et aLll2004l: 
iDekel fc Wo'3l2003L in SDSS and the Local Group), 
(m) AGN. Black hole m asses are correlated wit h their host 
spheroid properties (e.g. iTremaine et al.ll2002^ . The opti- 
cal AGN population, with high accretion rate and SFR, 
peaks near Ms,crit with little AGN activity at smaller 
masses, and is asso ciated with black- hole masses ^ 10* M© 
iKauffmann et al J 12003. in SDSS). Radio-loud AGNs, un- 
correlated with the optical activity and the SFR, dominate 
in larger haloes hosting ~ lO^^^M© black holes (G. Kauff- 
mann, private communication). 



1.2 Key Physical Processes 

The bi-modality imprinted on almost every global property 
of galaxies deserves a simple theoretical understanding. We 
propose that the main source of the bi-modality is the tran- 
sition from cold flows to virial shock-heating at a critical 
scale, in concert with feedback processes and gravitational 
clustering that emphasize the same characteristic scale. We 
address the cross-talk between these processes, and integrate 
them into a scenario which attempts to address simultane- 
ously the variety of observed phenomena. The key processes 
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(a) Cold infall vs. hot medium. The thermal behavior of 
the gas as it falls through the halo is qualitatively different 
below and above a critical mass scale of M ghock ~ 10^^ Mq 
jBirnboim fc DekeJl2003l: iKeres et 'ai]l2004l) . In less massive 
haloes, the disc is built by cold flows (~ 10''~^K), which 
are likely to generate early bursts of star formation. In more 
massive haloes, the infalling gas is first heated by stable 
shocks to near the virial temperature (~ lO^K). Near and 
above Mshock at 2 ^ 2 (and preferentially in isolated galax- 
ies) , streams of dense c old gas penetrate through the dilute 
shock-heated me dium jPardal et al.l I2001I : iKravtsovl 120031: 
iKeres et al.ll2004h (discussed in 33®. 

(b) Gravitational clustering. Non-linear gravitational 
clustering of the DM occurs on a characteristic mass scale, 
Mt , marking the typical haloes forming at a given epoch and 
the lower bound for groups of galaxies. We point out that 
the clustering scale, which varies rapidly with cosmological 
time, happens to coincide with Mahock at 2 ^ 1, and the 
interplay between these scales plays a role in determining 
the galaxy properties (0. 

(c) Feedback. We argue that the feedback processes affect- 
ing galaxy evolution are relatively ineffective near Mghock, 
largely due to the shock-heating process itself, and they 
therefore help emphasizing the imprint of this scale on the 
galaxy properties. Supernova and other feedback processes 
regulate star formation in the blue sequence below Mshock- 
Feedback by AGNs, or other sources, becomes efficient in 
haloes more massive than Mshock, because it preferentially 
affects the dilute shock-heated medium and may prevent it 
from ever cooling and forming stars, effects (JSJ- 

We show how the introduction of shock-stability physics 
crystallizes our understanding of the origin of the charac- 
teristic scales of galaxies. We argue that the combination 
of shock heating, feedback and clustering introduces a new 
feature in galaxy-formation modeling — a complete sup- 
pression of cold gas supply in haloes above a critical mass 
after a critical redshift. This could be the key to solving 
many of the open questions posed by the observations, fo- 
cusing on the bright-end truncation of the luminosity func- 
tion, the appearance of very red bright galaxies already at 
2; ~ 1 at the expense of big blue galaxies, and the indi- 
cations for massive starbursts at higher redshifts. We note 
that some of the issues have been addr essed in parallel, in 
a simi l ar spirit and in diffe r ent ways, by |BijnbQinifc^£ke 
|200,3[) . lBenson et al.l l|200,1) . lKeres et alJ (|2004 ) and lBinnev 
l|2004) . We make here a more thorough investigation into 
the cold flows versus shock heating phenomenon, relate it to 
the feedback and clustering processes, and attempt an inte- 
grated scenario that addresses simultaneously the variety of 
observed features. 

1.3 Outline 

In [21 we provide an improved present ation of the origi- 
nal analysis of spherical shock stability iBirnboim &: Dekell 
[2003 . hereafter BD03). In i|3|we compute the associated crit- 
ical mass scale in the cosmological context. In ®we describe 
how the phenomena is demonstrated in cosmological simula- 
tions, and learn about cold filaments in massive hot haloes 
at high redshift. In we address the cross-talk with the 
relevant feedback processes working alternatively below and 
above the critical scale. Then, in we integrate the above 



processes into a scenario which attempts to explain the ori- 
gin of the bi-modality and related features, and report first 
results from simulations that implement the new ingredi- 
ents. In ^we briefly discuss possible implications on other 
open issues in galaxy formation, and in we summarize 
our results, the proposed re-engineering of galaxy formation 
simulations, and the open theoretical issues. 



2 SPHERICAL SHOCK-STABILITY ANALYSIS 

The standard paradigm of disc formation 



jRees & Ostrikerl 


|1977| 




telM I1977I: IWhite & Reed Il978l: 


Bli^n 
I1991I 


lenthal. Faber. Prima9k & R.eesI Il984l: IWhite & FrenkI 
IMo. Mao & White!ll998l). which lies at the basis of 



all current models of galaxy formation, assumes that while 
a DM halo relaxes to a virial equilibrium, the gas that 
falls in within it is shock heated near the halo virial radius 
_Rv to the halo virial temperature. It is then assumed to 
cool radiatively from the inside out. As long as the cooling 
time is shorter than a certain global free-fall time (or the 
Hubble time), typically inside a current "cooling radius", 
the gas is assumed to accrete gradually onto a central disc 
and then form stars in a quiescent way. The maximum 
halo mass for efficient cooling was estimated to be on the 
order of ~ 10^^"^^Mq, and the common wisdom has been 
since then that this explains the upper bound for disc 
galaxies. However, early hints, based on one-dimensional 
simulations, indicated that this scenario cannot reproduce 
the sharp drop in the lum inosity function above this scale 
ijThoul fc Weinberd Even earlier studies, valid in 

the context of the pancake picture of galaxy formation, 
indicated tha t virial shock heating may not be as general 
as assumed (|Binnevl More advanced cosmological 

simula tions have started to reveal the presence of cold 
flows iFardal et alj I2OOII) . With the new data from big 
surveys such as SDSS, 2MASS and 2dF, and the detailed 
semi-analytic modeling (SAM) of galaxy formation, it is 
becoming clear that the observed scale is somewhat smaller 
and the drop is sharper than predicted by the original 
picture. It seems that the current models based on the 
standard paradigm have hard time trying to reproduce 
many of the observed bi-modality features summarized 
in This motivated us to attempt a closer look at the 
shock-heating mechanism. 

2.1 Spherical simulations 

Fig. shows the time evolution of the radii of Lagrangian 
gas shells in a spherical gravitating system consisting of gas 
(in this case with primordial composition) and dark matter, 
similar to the original simulations by BD03 using an accu- 
rate one- dimensional hydrodynamical code. Not shown are 
the dissipationless DM shells, which detach from the cos- 
mological expansion, collapse and then oscillate into virial 
equilibrium such that they deepen the potential well attract- 
ing the dissipating gas shells. The initial density perturba- 
tion was assumed to have a profile proportional to the linear 
two-p oint correlation function of ACDM (justified in iDekell 
[19811), and the final density profile roughly mimics the uni- 
versal NFW profile seen in cosmological simulations. The gas 
is cooling radiatively based on the atomic cooling function 
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3 4 5 

time [Gyr] 

Figure 1. Time evolution of the radii of Lagrangian gas shells 
(lines) in a spherical simulation of a protogalaxy consisting of pri- 
mordial gas {Z = 0) and dark matter. Temperature is marked by 
color. A shock shows up as a sharp break in the flow lines, namely 
a sudden slowdown of the infall, associated with an abrupt in- 
crease in the temperature. The lower discontinuity where the in- 
flow is brought to a final halt marks the "disc" radius, formed 
due to an artificial centrifugal force, (a) A massive system, where 
the virialized mass grows from 10" to 10^3 Mq. (b) A less mas- 
sive system, growing from 10^" to 10^^ Mq. A virial shock ex- 
ists only in systems more massive than a critical mass, while in 
smaller haloes the gas flows cold and unperturbed into the inner 
halo. With more realistic metallicities the critical mass becomes 



computed bv ISiitherland fc Dooi ta' ("1993"). The collapse of 
each gas shell is stopped at roughly 0.05-Rv by an artificial 
centrifugal force which mimics the formation of a central 
disc. 

The upper panel focuses on massive haloes of ~ 
10^^ Mq. As expected in the common picture, a strong shock 
exists near the virial radius, namely at roughly half the 
maximum-expansion radius of the corresponding shell. The 
virial shock gradually propagates outward, encompassing 
more mass in time. The hot post-shock gas is in a quasi- 
static equilibrium, pressure supported against gravitational 
collapse. The lower panel focuses on halo masses smaller by 
an order of magnitude, and shows an interesting new phe- 
nomenon. A stable shock forms and inflates from the disc 
toward the virial radius only after a total mass of more than 
a few times 10^^ M© has collapsed. In less massive systems, 
the cooling rate is faster than the compression rate required 
for restoring the pressure in the post-shock gas; had there 
been a shock, the post-shock gas would have become un- 
stable against radial gravitational contraction and unable 
to support the shock. In the specific case shown, with zero 
metallicity, the critical mass is biased low; with more re- 
alistic metallicities it becomes ~ 10^^ M© due to the more 
efficient cooling via metal lines (©. 



We demonstrate below (SjlJ that the behavior seen in 
the spherical simulations is robustly reproduced in cosmo- 
logical simulations. But first we wish to understand the 
spherical case, and use it for simple predictions. 



2.2 Post-shock stability criterion 

The existence or absence of a shock, as seen in the simula- 
tions, can be evaluated via a straightforward stability anal- 
ysis of the post-shock gas (first introduced in BD03). We 
provide here a brief, improved presentation of this analysis, 
followed by a more detailed estimate of the predicted critical 
mass for shock stability as a function of redshift (JJ^J. 

A stable extended shock can exist when the pressure in 
the post-shock gas is sufficient to balances the gravitational 
attraction toward the halo centre. The standard equation of 
state for an ideal gas expresses the pressure as P = (7 — 
l)pe where p and e are the gas density and internal energy 
per unit mass, and 7 — 5/3 for a mono-atomic gas. In the 
text-book case of no cooling, the adiabatic index is defined 
as 7 = (91n P/9 lnp)ad, and the system is known to be 
gravitationally stable once 7 > 4/3. When there is energy 
loss (e.g. by radiation) at a rate q per unit mass, we define 
a new quantity along the particle trajectories: 



d In P p q 

7off = -71 = 7 - - - 

amp p e 



(1) 



The second equality follows from energy conservation, e = 
—PV — q (where V = 1/p), plugged into the equation of 
state. Note that 7cff = 7 when g = 0. The difference be- 
tween the two is a ratio of characteristic rates for the two 
competing processes: the cooling, which reduces the pres- 
sure in the post-shock gas, and the compression due to the 
pattern of the post-shock infall, which tends to increase the 
pressure. If the compression rate is efficient compared to 
the cooling-loss rate, it restores the pressure necessary for 
supporting a stable extended shock, but otherwise the post- 
shock gas collapses inward under gravity, failing to support 
the extended shock. 

It is convenient to express the compression rate in the 
post-shock region as the inverse of a compression time, which 
we define by 



37 + 2 
7(37 - 4) 



21 



(2) 



with the factor V to be justified below, and the last equality 
referring to 7 = 5/3. For a spherical shock at radius rs, and 
a post-shock radial velocity Ui, we assume that the radial 
fiow pattern in the post-shock region is homologous, u/r = 
ui/rs. This is justified based on the spherical simulations 
described above, where the log-linear post-shock fiow lines 
in Fig. are nearly parallel straight lines. We then obtain 
using continuity 



tc 



Fr, 



(-V-u) (-3m) 



(3) 



The competing cooling rate in the post-shock region is 
expressed as the inverse of the standard radiative cooling 
time defined by 



f-cool 



(4) 
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where e = e{T) and q oc pA{T, Z), functions of temperature 
T and metallicity Z. Then in eq. Q 



7eff = 7 ■ 



tc 



(5) 



In order to test for stability, BD03 performed a pertur- 
bation analysis were the radius of a shell is perturbed by 
r — * r + (5r and the sign of the force, 5r/Sr, is computed. 
Writing Sr = u5t, using the homology, and assuming fur- 
ther that the gravity and pressure forces balance each other 
near the transition state between stability and instability, 
p~^VP = GM/r^, one obtains a restoring force, i.e. stabil- 
ity, for 

27 _ 10 
7 + 2/3 ^ 7 ■ 

The 7crit ~ 10/7 replaces the standard 7crit ~ 4/3 of the 
adiabatic case.^ 

Using eq. Q and the definitions of the time scales 
above, the shock stability criterion of eq. Q becomes the 
simple condition that the cooling rate should be slower than 
the compression rate: 



JcB > 7crit 



(6) 



> tc 



(7) 



Once the cooling rate is slower, the pressure gain by com- 
pression can balance the loss by radiative cooling, which 
allows the post-shock gas to be stable against global grav- 
itational collapse and thus support the shock. The factor 
r = 21/5 has been introduced in the definition of fcomp, 
eq. (|5J, in order to simplify this final expression. 

Note that the relevant quantity for stability is the ratio 
of rates associated with the two competing processes, inde- 
pendent of how slow each of them actually is in absolute 
terms. Each of the characteristic times could in principle be 
longer than the Hubble time - it is their ratio which deter- 
mines whether a stable shock is possible or the gas falls in 
subject to gravity, cold and unperturbed. 



2.3 Pre-shock quantities 

2.3.1 Compression rate 

Using the standard jump conditions across a strong shock, 
we can express the characteristic time scales (or 7eff) in 
terms of the pre-shock gas quantities. The jump condition 
for the radial velocity is 

7 + 1 

UO — Us = -{Ul — Us) , (8) 

7-1 

where Us is the radial shock velocity and uo is the radial 
velocity of the pre-shock gas. Then 



r(7 + l) 

3(7 - 1) \uo\ 



1 - 



(7- 



uo 



(9) 



^ If the spherical symmetry assumed above is replaced by planar 
symmetry, both for the shock and the gravitational field, the sta- 
bility criterion ■y^.g > 10/7 is replaced by 7eff > 10/11 (Birnboim, 
Dekel & Kravtsov, in preparation). One can therefore assume in 
general that the actual critical value lies somewhere between these 
two limits; if 7eff < 10/11 there is no stable shock, if 7(,fl > 10/7 
the conditions allow a stable shock, and if 10/11 < ■y^s < 10/7 
the shock stability depends on the local geometry. 



28 

5 \uo\ 



[l-3us 



5.48 Gyr--^(l-3us 



(10) 

(11) 



where its = Us/\uq\ (see H2.3.41 and the last expression as- 
sumes 7 = 5/3, Ts in 100 kpc, and Mo in 100 kms~^. If Us — 0, 
say, then tcomp is about 6 times larger than rs/|ito|, a typi- 
cal free-fall time from rs into the halo centre. At the virial 
radius, fcomp is comparable to the Hubble time at the corre- 
sponding epoch, but at inner radii it becomes significantly 
shorter. 



2.3.2 Cooling rate 

The cooling time (e.g. based on [Sutherland fc" Dopitalll993ll 



^cool — — 



(1 



(1 + ^ 



— pA(T, Z) 



(12) 



where A(r, Z) is the cooling function, k is the Boltzmann 
constant, e = nHe/nn, the mass per particle is m = prUp 
with jj, — {1 + 4e)/(2 + 3e), and the number of electrons per 
particle is x = (1 + 2e)/(2 + 3e). For 25% He in mass, one 
has e — 1/12, yielding n = 0.59. If we express the post- 
shock temperature as Te = T/IO^K, the post-shock baryon 
density as p-2s 
as A_22(T, Z) = A(T, Z)/10 



p/W ^^gcm ^, and the cooling function 



^^erg cm'' s ^ , we have 



tr. 



2.QlGyrpZlsTeAZl2{T,Z) 



(13) 



The post-shock gas density is related to the pre-shock 
density by the jump condition 



7+1 

7-1 



po = 4po , 



(14) 



and the post-shock temperature entering the cooling time is 
related to the pre-shock radial velocity uq via 

fcTi 2(7-1) , 2 3 2., ^-^2 

We note in passing that for a virial shock, where uq = 
— Vv (BD03), the post-shock temperature is actually 



Ti > - Tv 



(16) 



comparable to but somewhat smaller than the virial tem- 
perature as defined in eq. <A9ll . 



2.3.3 Stability criterion 

Using eas. llll and ll3l the critical stability condition becomes 

icool 
^comD 



P-lsTe A_L(T, Z) 

' 0.48 I Vw. o-\-i - 1 
rs\uo\ ^(1 — 3?is) ^ 



(17) 



with rs in 100 kpc and uo in 100 kms~^. Recall that eq. 1141 
relates p to po, and eq. 11511 relates T to mq. Thus, for given 
shock radius rs, shock velocity relative to infall Ug/|uo|, and 
pre-shock gas density po, once the metallicity Z is given and 
the cooling function A(T, Z) is known, one can solve eq. 1171 
for the critical values of T and the corresponding uq. When 
put in a cosmological context {^^, this solution is associated 
with a unique critical halo mass. 
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The stability criterion derived above, eq. © or eq. 
is found to work very well when compared to the results 
of the spherical simulations shown in Fig. Q When 7off (or 
tcooi/icomp) is computed using pre-shock quantities at a po- 
sition just outside the "disc" , we find that as long as the halo 
is less massive than a critical scale, before the shock forms, 
the value of 7cff is indeed well below 7crit and is gradually 
rising, reaching 7crit almost exactly when the shock starts 
propagating outward. The value of the y^s computed using 
the quantities just outside the shock then oscillates about 
7crit with a decreasing amplitude, following the oscillations 
in the shock radius seen in Fig. Q As the shock eventually 
settles at the virial radius, 7eff approaches 5/3, larger than 
7crit = 10/7, where the cooling is negligible. The same sta- 
bility criterion is found to be valid to a good approximation 
also in cosmological simulations i^^. 

2.3.4 Shock velocity 

What value of Ua is relevant for evaluating stability? In the 
inner halo, we use Us = 0. This is because, as the halo is 
growing in mass, the shock first forms in the inner halo and 
then propagates outward (Fig. 03). The onset of shock sta- 
bility is therefore marked by its ability to develop a velocity 
outward. 

During the stable phase when the shock is expand- 
ing with the virial radius, the spherical simulations indicate 
roughly ui ~ — Us (Fig.IIK), namely Us ^ 1/7 (eq.|HJ. This 
indicates that a small shock velocity of such a magnitude is 
appropriate for evaluating stability at the virial radius. 

Note that stability is harder to achieve when the shock 
is expanding relatively fast. In particular, in the extreme 
case Ua = 1/3, the post-shock velocity vanishes, iti = 
(eq.iJ. The compression rate becomes infinitely slow (eq.l^J, 
implying that the shock cannot be stabilized. 



3 SHOCK-HEATING SCALE IN COSMOLOGY 
3.1 Haloes in cosmology 

We wish to translate the critical stability condition, eq. ((7|l 
or eq. 11711 . into a critical post-shock temperature, and the 
corresponding critical halo virial velocity and mass as a func- 
tion of redshift. Eq. 11711 has a unique solution when com- 
bined with the two virial relations between halo mass, veloc- 
ity and radius (eg. lAGL and the relation between post-shock 
temperature and pre-shock infall velocity (eg. llSH . 

As summarized in Appendix A, the time dependence of 
the virial relations can be expressed in terms of the conve- 
nient parameter 

^ = (AaooS^mo.afto.T)"'''^' 



(18) 

where a = 1/(1 + z) is the cosmological expansion factor 
and the other parameters are of order unity. The parame- 
ters f2mo.3 and ho.7 correspond to today's values of the cos- 
mological mass density parameter and the Hubble expansion 
parameter respectively, and for the standard ACDM cosmol- 
ogy adopted here they are both equal to unity. The param- 
eter A200 is the virial density factor given approximately in 
eq. lA8t : at redshifts 2: > 1 it is A200 — 1, but at lower 
redshifts it becomes somewhat larger, reaching A200 — 1-7 
at z — 0. 



3.2 Compression rate 

For a shock at the virial radius, Vs = Rv, we approximate 
uo = —Vv, as predicted by the spherical collapse model in 
an Einstein-de Sitter cosmology (BD03, Appendix B). 

When the shock is at an arbitrary inner radius r, where 
the infall velocity is \u\, we multiply Rv and Vv by appropri- 
ate factors fr = r/Rv and /u = |u|/Vv (discussed in H3.5II . 
Then eq. lllll becomes 



tc 



U.3GyrA'/\frf-\l-3uay 



(19) 



3.3 Cooling rate: gas density 

In order to express the cooling time of eq. 1131 in terms 
of cosmological quantities, we first evaluate the pre-shock 
baryon density, which we write as 



Pb = 4 /b (p/p)vir A pu fp 



(20) 



Here pu is the universal mean mass density (eg. IA4II . and 
A is the top-hat mean overdensity inside the virial radius 
feg. lASt . The factor (p/p)vir translates p, the mean total 
density interior to J?v, to p, the local total density at R^,. 
The effective baryonic fraction /b turns it into a pre-shock 
baryonic density. The factor 4 stands for the ratio between 
the post-shock gas density and the pre-shock gas density.^ 
The factor fp = p{r)/p{Rv) reflects the ratio of the actual 
gas density at some radius r within the halo to its value at 
the virial radius (see H3.5t . 

The ratio p/p at the virial radius is derived for the 
universal NFW halo density profile revealed b y cosmologi- 
cal simulations iNavarro. Frenk fc Whit^ll997^ . For a virial 
concentration parameter c, this ratio is 



3(1 -^c)2 



ln(l + c) 



(1 + c) 



(21) 



A typical concentration of c = 12 is associated with p/p — 
0.17; we therefore express the approximate results below us- 
ing the factor fp,o.i7 = (p/p)/0.17. In our more accurate 
evaluation of the critical scale ( H3.7L we model the depen- 
dence of the mean concentration on mass and time using the 
fit of Bullock ct al. ( 200J() for the ACDM cosmology: 



c(M,a) = 18Mri°-"a 



(22) 



The effective baryon fraction /b may in principle be 
as large as the universal fraction ~ 0.13, but it is likely 
to be smaller because of gas loss due to outflows. For the 
approximate expressions we define /b.o.os = /b/0.05. 

The gas density at r, eq. 12011 . thus becomes 



P-28 =0.190A-Vb,.05/p/p,0.1 



(23) 



Inserting this baryon density into eq. 11311 . the cooling time 
becomes 



tc 



13.7 Gyr AVr.os/p" VpAiT ^ AZ^T, Z). 



(24) 



In the spherical simulations, the relevant factor relating the 
baryon density to the DM density in eq. l2Ui is actually closer to 
~ 3 because of a "bump" in the DM density just inside the virial 
radius. 
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The coolinK function that we use below (based on 
ISutherland fc Dopit a 1993) neglects two physical processes: 
Compton scattering off the cosmic microwave background 
and the possible effect of external radiation on the cool- 
ing rate through the reionization of Hydrogen. Based on 
the more complete cooling function as implemented by 
iKravtsov fc GnedinI i2004fl . one learns that these processes 
become important only for densities below ~ 10"'^* and 
~ 10~^® gcm~^ at z ~ and 4 respectively. Using eq. i'2'6l . 
we conclude that while these processes may have a certain 
effect on the cooling rate near the virial radius, they should 
be negligible once the analysis is applied inside the inner 
half of the halo, where the critical scale for shock heating 
is determined in practice. We address these effects in more 
detail elsewhere (Birnboim, Dekel fc Loeb, in preparation). 



3.4 Metallicity 

The metallicity near the virial radius and in the inner halo, 
which also enters the cooling rate, is one of our most uncer- 
tain inputs. For the mean metallicity Z (in solar units) as a 
function of redshift z we use the two-parameter functional 
form 



logiZ/Zo) 



(25) 



where Zq is today's metallicity and the slope s governs the 
rate of growth. 

An upper limit may be imposed by the hot. X-ray emit- 
ting Intra-Cluster Medium (ICM) at low redshifts, which 
indicate Zq ~ 0.2 — 0.3. The ICM metallicity evolution in 
Semi- Analytic Models, assuming a range of different feed- 
back recipes, yields consistently an average enrichment rate 
of roughly s ~ 0.17 (R. Somerville, p rivate communication; 
iDe Lucia. Kauffmann fc Whitel2o'Ml . We adopt this enrich- 
ment rate s in our modeling below. 

A realistic estimate of the metallicity near the virial ra- 
dius (or perhaps a lower limit for the inner halo) may be 
provided by Civ a bsorbers in the Int er-Galactic Medium 
(IGM) at z ~ 2 - 4 iSchave et al.ll2003ft . At densities appro- 
priate to typical NFW haloes at z = 3 (with c = 3), namely 
Pvir — 53pu, they measure an average of [C/H] = —2.47. 
Silicon measurements, on the other hand, seem to indicate 
a metallicity that is about five times larger (A. Aguirre, 
private communication). If one takes the geometrical mean 
between the metallicities indicated by Civ and by Si one 
has Z{z = 3) ~ 0.0075. This translates to Zo = 0.025 ff 
s = 0.17. 

We note that another popular indicator, Mgn, indi- 
cates consistently Z < 0.01 within 50 — 100 kpc of galaxy 
cen tres at z ^ 1 (pri vate communication with J. Charlton; 
e.g. iDing et alJboOSD . 

The Damped Lyman-Alpha Systems (DLAS) are be- 
lieved to sample cold gas deeper inside the haloes, and can 
thus provide another interesting limit. Observations in the 
range z = 1 — 4 Iprochaska ct al. 2003) can be fitted on aver- 
age by eq. I25II with Zq ~ 0.2 and a somewhat steeper slope 
s ~ 0.26. However, a fit with s = 0.17 (and then Zq = 0.1) 
is not ruled out. 

Based on the above estimates, we adopt as our fiducial 
metallicities Zo — 0.03 at i?v and Zo =0.1 at the "disc" 
radius ~ 0.1-Rv, both with an enrichment rate s = 0.17. 



3.5 Inside the halo 

For a shock in the inner halo we wish to estimate the factors 
fr, fu and fp. 

Empirically from the spherical simulation of BD03, for 
a shell encompassing a mass just shy of the critical mass (as 
well as from the toy model of BD03 of gas contracting in a 
static isothermal sphere), we estimate for fr = r/Rv 



f, 



Po{r) 



J r ! 
r-2.1 
J r 7 



r ~ O.liiv 



(26) 



po(-Rv) 

We adopt below fp = f^^ at /,. = 0.1. 

Energy conservation assuming pure radial motion inside 
a static singular isothermal sphere yields 

uo{r) 



= [l + 2/b(/-^-l) + 2(l-/b)ln/. 



-111/2 



(27) 



Uo(i?v) 

For fr — 0.1 and /b ~ 0.05 this gives the estimate f^ ~ 2.5. 

Based on the definition of /„, the temperature behind 
a virial shock is related to the temperature obtained from 
the stability condition at radius r by 



T{R.) = f-' T{r) 



3.6 Crude explicit estimates 



(28) 



The critical temperature for stability, as obtained by com- 
paring tcooi and tcomp in the cosmological context, eq. 1241 
and eq. 11911 . is 



Te\Zl2iT,Z) = 1.04^^^/^^ , 
where 

F = frfn^fp /b,.05/p,0.17(l - 3fis)"^ 



(29) 
(30) 



The cooling function as computed by 
ISutherland fc Dopital (|l993h can be crudely approximated 
in the range 0.1 < Te < 10 by 



A_22 ^ 0.12 Zg,o^3 TfT' + 0.02 



1/2 



(31) 



where Z0.03 = ^/0.03, with Z in solar units. The above ex- 
pression is valid for — 2.5 ^ log Z ^ 0, and at lower metallici- 
ties the value of A is practically the same as for log Z = —2.5. 
This fit is good near Te ~ 1 for all values of Z. The first 
term refers to atomic cooling, while the second term is due 
to Bremsstrahlung. For an approximation relevant in haloes 
near Afshock we ignore the Bremsstrahlung term, which be- 
comes noticeable only at higher temperatures. One can then 
obtain in eq. 12911 an analytic estimate for the critical tem- 
perature: 



Te -0.36^1-=*/'' (4.03i^)'^' 



(32) 



Using eq. 11511 and eq. 128L with jitoj = Vv just outside 
the virial radius, we then obtain for the critical velocity and 
mass 

V,oo ^ 1.62^-^/« (Z^dsFY^^ f-' (1 + u.)-' , (33) 

Mn ^ 25.9A^/' [Z^^Ff^^ f'' (1 + u,)-'' . (34) 

A comment regarding the Us dependence of our results. 
The critical temperature depends on the shock velocity Ma 
via _F, T oc (1 — 3{ts)~^^^, refiecting the Us dependence of 
tcomp. The critical temperature is thus monotonically in- 
creasing with Us- An additional Us dependence enters when 
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the temperature is translated to a critical virial velocity us- 
ing the jump condition, V (x {1 + Ua)~^T^^^ , and then to 
a critical mass, Af oc (1 + fia)"^T^/^. For a slowly moving 
shock, Ws <C 1/3, the combined Ma dependence of the critical 
mass is M oc [l+(9/4)usl(l— Sus) ^ 1— (3/4)ms — a decreas- 
ing function of Us- This means that at a given radius in a 
halo of a given mass, when everything else is equal, a slowly 
expanding shock is actually more stable than a shock at rest. 
For example, if the shock is expanding with = 1/7 rather 
than Us = 0, the critical mass is smaller by about 24%. How- 
ever, recall that stability is harder to achieve when the shock 
is expanding relatively fast, and the compression completely 
vanishes if ^ 1/3 ( 12.3.41 . 

For actual crude estimates of the critical scales at z — 0, 
we assume /b,o.05 — fpfi.n — 1. For a shock at the virial 



10' 



radius, fr = fu 
1/7, and obtain 

n ^ 0.5, Vioo 



fp = 1, we assume Zq — 0.03 and Ub 



1.6, Mil ^ 26. 



(35) 



At an inner radius closer to the disc vicinity, say fr = 
0.1, we estimate /„ ~ 2.5 and fp ~ 100 (> l3.5t . Assuming 
Zo ~ 0.1 and Us — 0, we obtain 



Te ~ 1.1, Vio 



1.1, Afii 



8.8. 



(36) 



We see that for a shock at r ~ 0.1-Rv, the expected critical 
mass is smaller than at i?v, somewhat below ~ W^^Mq. 

The above estimates are useful for exploring the quali- 
tative dependences of the critical values on redshift, metal- 
licity and gas fraction. For example, in eq. (1341 . the explicit 
redshift dependence and the decrease of metallicity with red- 
shift tend to lower the critical mass toward higher z. On the 
other hand, the decrease of halo concentration with z (i.e. 
increase of /p.o.ir), and the possible increase of the effective 
gas fraction with z (®, tend to push the critical mass up 
at higher z. 

3.7 More accurate estimates 

We now obtain a better estimate of the critical tempera- 
ture (and then critical mass and virial velocity) by solving 
eq. 1291 using the exact cooling function of SD93 and adopt- 
ing specific models for the time evolution of metallicity and 
halo structure. The results are presented graphically. 

The baryon density is computed assuming an NFW pro- 
file whose concentration evolves in time as in eq. (1221 . The 
effective fraction of cold gas is assumed to be /b = 0.05, 
motiv ated by best fits of semi-analytic models to the Milky 
Way (|Klvpin. Zhao fc Somervilld l2002^ and by fitting the 

ACDM halo mass function to the observed luminosity func- 

I 11 TJs 

tion near lIBell et al . 2003). The metallicity evolution is 

parametrized as in eq. (1251 with s = 0.17 for today's metal- 
licities in the range Zq = 0.03 — 0.3. Upper and lower esti- 
mates for the critical scales are evaluated at the virial radius 
and at r = 0.1-Rv respectively, using the crude estimates of 
H3.5I In the following figures the shock is assumed to be at 
rest. Us = 0. 

Figure|21shows the critical mass as a function of redshift. 
At a typical inner-halo radius, r = 0.1 -Rv, we plot the curves 
for three different current metallicities: Zq = 0.03, 0.1, 0.3. 
The critical halo mass, for Zo = 0.1, is ~ 6 X 10"Mo 
quite independent of redshift. The uncertain metallicity in- 
troduces a scatter by a factor of 2 up and down (for z < 2.5). 
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Figure 2. Critical shock- heating halo mass as a function of 
redshift. The three solid (red) curves refer to a shock at the inner 
halo, r = 0.1-Rv, with different metallicities as indicated; the 
middle curve (Zq = 0.1) is our best estimate. The dash-dotted 
(magenta) curve refers to a shock at the virial radius with Zq = 
0.03. The other parameters used are: /b = 0.05, Us = 0., s = 0.17 
(see text) . Shown for comparison (short dash, blue) are the Press- 
Schechter estimates of the forming halo masses, corresponding 
to 1-(T (M*) and 2-a, where the fractions of total mass in more 
massive haloes are 22% and 4.7% respectively. Also shown is the 
critical mass for supernova feedback discussed in J5] (long dash, 
green) . 
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Figure 3. Same as Fig. l2lbut for the corresponding halo virial 
velocity. 



An upper limit of ~ 2 x 10 M0 is obtained for a shock 
at _Rv when a correspondingly low metallicity is assumed, 
Zo = 0.03. When the assumption of its ~ is replaced by 
■its — 1/7, allowing the shock to expand with the virial radius 
as seen in Fig. the critical mass at 7?v with Zo ~ 0.03 
becomes comparable to that at 0.1 with Zo = 0.3. 

Figure |3 shows the corresponding virial velocity. At 
z = the critical virial velocity for a shock in the inner halo 
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metallicity Z [Zq] 
Figure 4. Critical siioclc-lieating halo mass as a function of 
metallicity at redshift 2 = 0. The solid (red) curve refers to a 
shock at the inner halo, r = 0.1 -Rv The dash-dotted (magenta) 
curve refers to a shock at the virial radius. 



the corresponding fractions of the total mass encompassed 
in haloes exceeding the mass M are 22% and 4.7% respec- 
tively. One can see in Fig. |21 that Mshock coincides with M* 
at z ~ 1, and with the 2-a mass at z ~ 3.4. By z ^ 2, 
say, most of the forming haloes are significantly less mas- 
sive than Mshock. When embedded in a large-scale high-cr 
density peak, the distribution of forming haloes at a given z 
may shift toward more massive haloes. In fact, the most mas- 
sive halo in a volume of comoving size ~ lOOMpc is likely 
to be more massive than IO'^^Mq at all relevant redshifts 
[z < 6, say). Nevertheless, the qualitative result concerning 
the majority of the haloes remains valid. We conclude that 
m the vast majority of forming galaxies the gas has never 
been shock-heated to the virial temperature - it rather flows 
cold all the way to the disc vicinity. 

We note that the values obtained for Mshock at 
low redshifts are compatible with the observed bi- 
modality/transition scale summarized in ^ The estimates 
in the inner halo, using the lower and upper limits for 
Zo, indeed border the observed characteristic halo mass of 
~ 6 X IO^^Mq. The upper-limit estimate at _Rv corresponds 
to a halo mass similar to that of the Milky Way. 



is ~ 120 km s^ , with a ±30kms^ scatter due to metallic- 
ity. The critical virial velocity increases monotonically with 
redshift, to ~ 200kms~^ near 2 ~ 3 (a crude fit to the 
redshift dependence is Vv ~ 120 + 28z). 

The dependence on metallicity at z = is highlighted 
in Fig. 01 The metallicity enters strongly through the cool- 
ing function A(T,Z). The critical mass grows roughly like 
Z^^^, as predicted in eq. 13411 . so it spans about an order of 
magnitude over the whole metallicity range. 

The analytic estimates of eq. 1331 and eq. 13411 . based 
on the approximate cooling function, turn out to provide 
good estimates in most cases, and can therefore be used for 
extending the results analytically to any desired choice of 
the relevant parameters. 

We learn that the critical halo mass for shock stabil- 
ity at the disc vicinity, Mshock('"disc), is somewhat smaller 
than for a shock at the virial radius, Mghock(-Rv). This re- 
sult is robust: it is true even if the metallicity at the virial 
radius is smaller by an order of magnitude than the metal- 
licity at the disc, and even when Us at _Rv is as large 
as 1/7. This means that as the halo is growing in mass, 
the conditions for a stable shock develop first in the in- 
ner halo and somewhat later in the outer halo. Thus, in 
haloes of mass M < Mshock (j'disc), we expect cold flows 
with no shock heating throughout the halo. In the other ex- 
treme of haloes of mass M > Mshock (J?v), we expect shock 
heating of most of the gas by a shock near the virial ra- 
dius. In haloes of mass in the narrow intermediate range 
Mahock(rdisc) < M < Mshock(-Rv), we expect shock heating 
somewhere inside the halo, preventing most of the gas from 
falling in and giving rise to a hot medium. This predicted 
mass range, of a factor of 2-3 in mass, is consistent with 
the range of transition from all cold to mostly hot seen in 
cosmological simulations O. 

Also shown in Fig. |5| are the typical masses of haloes 
forming at different redshifts, the 1-a (termed M») and 2- 
a halo masses according to the Press- Schechter formalism, 
eq. lAlSt . According to the improved Sheth-Tormen version, 



4 COLD STREAMS IN HOT HALOES 
4.1 Cosmological simulations 

Cosmological hydro simulations indicate that the phe- 
nomenon of cold flows is a general phenomenon not re- 
stricted to spherical symmetry. Fig. |3 displays snapshots of 
an Eulerian simulation from Birnboim, Dekel, Kravtsov & 
Zinger (in preparation, BDKZ05).* Shown are maps of gas 
temperature in two epochs in the evolution of a protogalaxy: 
one at 2 ~ 4, when the halo is already relative massive, and 
the other at 2 ~ 9, when the halo is still rather small. While 
the more massive halo, near the critical scale, shows a hot 
gas component near the virial temperature behind a virial 
shock, the smaller halo shows only cold flows inside the virial 
radius. 

Similar res ults have been obta ined earlier from SPH 
simulations by iFardal et al.l i200ll) . who emphasized the 
feeding of galaxies by cold flows preferentially at early 
epochs. Based on our spherical analysis, we understand that 
this redshift dependence mos tly reflects t he smaller masses 
of haloes at higher redshifts. iKeres et alJ (|2004) have ana- 
lyzed similar SPH simulations and presented the case for the 
two modes of infall, cold and hot, in more detail.^ For exam- 
ple, in their figures 1 and 2 they demonstrate the two-mode 
phenomenon very convincingly by showing the distribution 
of particles and their trajectories in temperature-density di- 
agrams. Their most informative Fig. 6 shows the fractions 
of cold and hot infall as a function of halo mass at different 
redshifts. For all haloes of masses below a critical mass the 
infall is predominantly cold. Near the critical mass there is a 
relatively sharp transition into a hot mode, which becomes 

^ A d escription of the simulation technique can be found in 
iKravts ov (2003), where it was used for other purposes. 
^ The hot phase becomes an "infall" mode in this simulation after 
the gas cools, but in reality it may be kept hot and be prevented 
from falling in by feedback effects, J5] 
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Figure 5. Snapshots from a cosmological hydrodynamical simu- 
lation (Birnboim, Dekel, Kravtsov & Zinger, in preparation, here- 
after BDKZ05) showing the gas temperature in a slice of a proto- 
galaxy at two different epochs, when it has two different masses. 
The temperature is relative to the virial temperature of the halo 
at that time. The side of each slice is scaled to be 3-Rv (numbers 
in comoving h~^Mpc). Top: At z ~ 4, when the halo is already 
relatively massive, M ~ 3x 10^^ Mq. Much of the gas is heated by 
a strong shock near the virial radius (white circle). Cold streams 
penetrate through the hot medium deep into the halo. Bottom: 
At z ~ 9, when the halo is still small, M ~2x 10^° Mq. The gas 
flows in cold (T <C Tv), showing no evidence for shock heating 
inside the virial radius (circle). 



dominant above the critical mass. The transition from 100% 
cold to more then 50% hot occurs across a range of only a 
factor ~ 2 in halo mass. In this simulation, where zero metal- 
licity is assumed, the transition mass is M ~ 3 X 10"Mo 
at all redshifts in the range — 3. A similar transition as 
a function of halo mass, and the constancy of the critical 
mass as a function of redshift, are both reproduced in the 
Eulerian cosmological simulations studies in BDKZ05. 

The spherical simulations described above ( H2.1II . and 
the corresponding analytic analysis (321 yield very sim- 
ilar results. In fact, our analytic predictions for the case 
of zero metallicity match the critical mass measured by 
iKeres et al.l i2004ll remarkably well. 

The stability criterion derived in the spherical case 
turns out to be valid locally in the cosmological simulations 
where the non-spherical features are pronounced. BDKZ05 
use this criterion to identify the cold streams and hot me- 
dia in the simulations without explicit information concern- 
ing the presence or absence of actual shocks. When test- 
ing the criterion in these simulations, in which the hot and 




Figure 6. A snapshot from a cosmological hydrodynamical sim- 
ulation (Zinger, Birnboim, Dekel & Kravtsov, in preparation, 
ZBDK05) showing the gas density of the cold flows (T < 10'''^K, 
p > 10~^^''^ g cm~^) within the virial radius of the same massive 
galaxy shown in Fig.|F]at z = 4. The cold phase is filamentary. In 
the outer radii, the gas filaments tend to ride on the large-scale 
DM filaments feeding the halo. About half the mass of the cold 
phase is in dense clumps. 



cold phases may be present in the same halo, the local gas 
properties at each position is first transformed to post-shock 
quantities, as if there was a shock there, and the stability 
is evaluated based on the derived value of jcti there. The 
resultant maps of "f^g resemble quite well the temperature 
maps of the actual simulation. This demonstrates that the 
wisdom gained by the spherical analysis is applicable in the 
general case. 



4.2 Filaments in the simulations 

The cosmological simulation of Fig. |S] also reveal that mas- 
sive haloes at high redshift can have cold streams embedded 
in the hot medium. While the hot medium is rather spher- 
ical, the cold fiows are filamentary and sometimes clumpy. 
Fig. |S] highlights the filamentary nature of the cold gas in 
the same b ig halo. This pheno menon is consistent with the 
findings of iKeres et al.l i2004) . that the cold mode may in 
some cases co-exist with the hot mode even above Afghock, 
especially at z > 2 (their Fig. 6) and preferentially in rel- 
atively isolated galaxies. We wish to understand the origin 
of this phenomenon, and learn about its dependence on cos- 
mological time. 

' ■ ■ simulation results of iKeres et alJ i2004l) provide 
several additional relevant clues. First, their figures 17 and 
18 indicate that the cold infall indeed tends to be filamen- 
tary, especially at high redshifts, while the hot mode is more 
spherical. They report that the directional signal measuring 
filamentary infall in the cold accretion mode is stronger for 
haloes above Afshock while the infall is more isotropic below 
it. 

Second, iKeres et alJ i2004l display in their figure 13 
the environment dependence of the gas infall modes, show- 
ing that the cold and hot modes dominate at low and high 
neighborhood densities respectively. We learn from this fig- 
ure that at z < 2 the cold mode dominates for galaxy densi- 
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ties below rigai ~ 0.3( ft~^Mpc)~^ and becomes negligible at 
larger environment densities, while at 2 = 3 the cold mode 
is more pronounced than the hot mode for all neighborhood 
densities up to rigai ~ 10( /i~^Mpc)~^. The correlation be- 
tween the environment density and the host halo mass im- 
plies that this z variation of the environment dependence 
could be partly attributed to the finding of a significant 
residual cold mode in massive haloes at high z (e.g., their 
figure 6). 

Third, iKeres et af] (|2004|) show in their figure 16 that 
the cold accretion is on average of higher density than the 
hot mode. This is by only a factor ~ 2 (probably underesti- 
mated because they mix small and large haloes), but since 
the shock is responsible for a density increase by a factor of 
4, the actual overdensity of the hypothet ical post-shock (or 
pre-sh ock) gas is more like ~ 8. Similarlv. lNagai fc Kravtsovl 
i2003f) find in their simulation of a massive halo that the fila- 
mentary structure is associated with gas entropy (oc T / p^^^) 
far below that of the surrounding halo gas. In the high- 
resolution simulation shown in Fig. ^ we find that the den- 
sity in the cold streams is actually higher than the surround- 
ing gas density by two orders of magnitude or more. The 
higher gas density in the filaments is associated with a more 
efficient cooling which prevents a shock from forming along 
the filaments (^. 

We find that the cold gas filaments at the halo outskirts 
are strongly correlated with the DM filaments that feed this 
halo (reported in detail in Seleson & Dekel, in preparation, 
hereafter SD05). These filaments are part of the large-scale 
cosmic web. They enter massive halos at high redshift as 
narrow streams with a density higher than the halo average 
by a factor of a few. The initial overdensity of the gas flowing 
along the DM filaments scales with the DM density, while 
its infiow velocity is comparable to the halo virial velocity. 
As a result, the initial cooling time in the thin filaments is 
shorter by a factor of a few than in the surrounding spherical 
halo, while the compression time is comparable in the fila- 
ments and the host halo. Eq. O then implies that the thin 
filaments have a harder time supporting a stable shock. The 
gas filaments remain cold, and become denser as the stream 
penetrates through the hot medium into the halo center. The 
result is that in massive haloes at high redshift the critical 
halo mass for shock heating in the filaments feeding them is 
larger than the estimate for a spherical virial shock derived 
in ^ We provide below a crude estimate for this revised 
critical mass. 



4.3 Interplay with the Clustering Scale 

What is the reason for the appearance of cold streams in 
massive haloes above Afghock at high z? First recall that 
there is another characteristic scale in the problem — the 
scale of nonlinear clustering M*, determined by the shape 
and amplitude of the initial fluctuation power spectrum and 
its growth rate. The masses for 1-a haloes (M*) and 2-a 
haloes, based on eq. lAlSt with v = 1 and 2 respectively, are 
shown again in Fig.Q We note that Mshock ~ M, at z ^ 1, 
while A^ahock ^ at z > 2. This means that ~ lO^^A'/0 
haloes are typical at z < 1 but they are the highest rare 
peaks at z > 2. We argue that this is responsible for the dif- 
ference in the cold-filament behavior of Af ^ Afshock haloes 
in the two epochs. Since the large-scale structure of dark 
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Figure 7. Cold streams and shock-heated medium as a function 
of halo mass and redshift. The nearly horizontal curve is the typ- 
ical threshold mass for a stable shock in the spherical infall from 
Fig. El below which the flows are predominantly cold and above 
which a shock-heated medium is present. The inclined solid curve 
is the upper limit for cold streams from eq. 1401 with / = 3; this 
upper limit is valid at redshifts higher than z^rit ~ 1 — 2, defined 
by Afgijock > fMt. The hot medium in M > Mohock haloes at 
> 2crit hosts cold streams which allow disc growth and star 
formation, while haloes of a similar mass at 2 < z^rit are all hot, 
shutting off gas supply and star formation. 



matter is roughly self-similar in time (when measured in 
terms of Af* and the background universal density), we can 
learn about the difference between typical and rare haloes by 
comparing M ~ Af* and M 3> Af* haloes in a single simula- 
tion snapshot. One can see in any high-resolution cosmolog- 
ical N-body simulation (e.g., the "Millennium Run", visu- 
alized in www.mpa-garching.mpg.de/galform/millennium) 
that the rare massive haloes tend to be nodes fed by a few in- 
tersecting relatively narrow filaments which are denser than 
the virial density of these haloes. On the other hand, a typ- 
ical Aft halo is commonly embedded in a single filament of 
the cosmic web, and this halo thus sees a wide-angle in- 
fiow pattern in which the matter density is comparable to 
the virial density (this is quantified in SD05). This explains 
why Si lO'^^AfQ haloes are fed by narrow dense streams at 
z > 2 but not at z < 1. 

A crude way to estimate the maximum halo mass for 
cold streams at a given redshift is as follows. Recall that the 
critical ratio for shock stability f eg. 1171 scales roughly like 



^cool P 

oc 

^comp 



(37) 



with T, R and V the halo virial quantities and p the gas 
density. At a given epoch, the typical halo p and R/V are 
roughly independent of halo mass (based on the definition of 
the virial radius), so with the virial relation T oc M^^^ , and 
approximating the cooling function with A oc feo. I31II . 
the critical ratio for spherical infall in virialized haloes is 



'=°""P / halo 



/ Af \-^/^ 

V Afshock / 



(38) 
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The analogous critical ratio in the dense streams inside 
a halo of mass M, assuming that RV~^ in the streams is the 
same as in the halo, is inversely proportional to the density 
enhancement Pstroam/phaio (eg. I37II . Our estimates from N- 
body simulations indicate that Pstrcam/pvir ^ (/M*/M)~^'^^ 
with / ~ 3 (SD05). With eq. I^Hl one obtains 




For this ratio to equal unity in the streams, the critical halo 
mass is 



Met ream /, r Mshock , /M. < Mshock . (40) 

This maximum mass for cold streams is shown in Figure Q 
At low z, where /M* > Mshock, cold streams exist only for 
M < Mshock- At high z, where /M, < Mshock, cold streams 
appear even in M > Mshock haloes where shocks heat part 
of the gas, as long as M < Mstroam- The critical redshift z^it 
separating these two regimes is defined by 

/M.(Zcrit) = Mshock . (41) 

This scenario is consistent with the cosmological hy- 
drodynamical simulations. The shock-heating mass explains 
the transition from cold to hot at a given mass roughly inde- 
pendent of z, and the presence of cold streams above Mshock 
at « > Zcrit explains the dependence of the cold mode on 
redshift and environment. Besides its dependence on halo 
mass, the environment effect (e.g., iKeres et alJl20o3l may 
also be due to the survivability of cold streams in different 
environments. While streams could survive unperturbed in 
relatively isolated galaxies, they are likely to be harassed by 
the active intergalactic environment in dense groups. The 
environment dependence may therefore also reflect varia- 
tions in the HOD at a given halo mass. The properties of 
cold flows in haloes as a function of halo mass, redshift, and 
grouping deserve a detailed analysis using high-resolution 
cosmological hydro simulations. 



5 FEEDBACK & LONG-TERM SHUTDOWN 

Once the halo gas is shock heated in massive haloes, what 
is the process that keeps it hot and maintains the shutdown 
required by the bi-modality? Is it also responsible for the 
rise of M/L with mass above Ms.crit (and the absence of 
cooling flows in clusters)? Several feedback mechanisms can 
heat the gas. We suggest that they have a minimum effec- 
tiveness in haloes ~ A/shock. This can be largely induced 
by the shock heating itself, and in turn it can amplify the 
bi-modality features. Some of the feedback mechanisms are 
limited to smaller haloes, while others, such as AGN feed- 
back, are likely to be important in more massive haloes. The 
latter can be triggered by the shock heating and then help 
maitaining the gas hot for a long time. 

5.1 Below the shock- heating scale 

(a) Supernova feedback. Based on the physics of super- 
nova (SN) remnants, the energy fed to the gas in haloes of 
r ~ lO^K is proporti onal to the stellar m ass despite signif- 
icant radiative losses jDekel fc Silklll98d) . When compared 



to the energy required for significantly heating the gas, one 
obtains a maximum halo virial velocity for SN feedback, 
VsN — 120kms~^. This is only we akly dependent on the 
gas fraction, density or metallicity jDekel fc SilklllQsH eq. 
49), and is therefore insensitive to redshift. Only in potential 
wells shallower than Vsn can the SN feedback significantly 
suppress further star formation and regulate the process. 
Figure |3 shows Vsn and Figure |5| shows the correspond- 
ing mass versus redshift. With an effective /b ^ 0.05, the 
corresponding stellar mass at 2 = is ~ 3.5 X 10^° Mq, 
practically coinciding with the bi-modality scale. The sim- 
ilarity of the SN and shock-heating scales is partly a coin- 
cidence, because the nuclear origin of the initial SN energy 
has little to do with galactic cooling or dynamics. Ifowever, 
there is an obvious similarity in the cooling processes and 
in the asymptotic behavior of a SN remnant, which is not 
a strong function of its initial energy. The distinct correla- 
tions between the properties of galaxies below Ms,crit indeed 
point at SN feedback as its primary driver. These correla- 
tions define a "fundamental line", V oc M°''^, Z oc M°''*, 
H PC M°'^, where n is surf ace brightness jKauffmann et aD 
l2003l:lTremonti et all2004^ SN feedback can ex plain the ori- 
gin of the fundamental line iDekel fc Wool2003^ based on (a) 
the above energy criterion, which implies Ms/M oc V^, (b) 
the virial relations (ea. lA6L (c) the instantaneous recycling 
approximation, Z oc Ms/Afgas, and (d) angular- momentum 
conservation, Rt o c A_R, with A a constant spin parameter 
jFall fc Efstathioul UQSOl . 

(b) UV-on-dust feedback. Also working below Afshock 
are Momentum-driven winds due to radiation pressure on 
dust grains, arising from the continuum absorption and 
scattering of UV photons emitted by starbursts or AGNs 
iMurrav. Quataert fc ThompsonI (|20^). The dust survives 
and can provide sufficient optical depth if the gas is cold and 
dense, e.g., in the cold flows below Afshock, which can also 
provide the starbursts responsible for the required UV flux 
and metals. Since dust grains cannot survive above ~ lO^K, 
Aishock imposes an upper bound for this feedback. 

(c) Photoionization feedback. The UV f rom stars and 
AGNs ionizes most of the g as after z ~ 10 llBuUock et alJ 
l2000l:lLoeb fc Barkanall200lli . heats it to > lO'^K, and pre- 
vents it from falling into haloes below the Jeans scale Vv ~ 
30kms^^ (Gnodin 2000). As the ionization persists for cos- 
mological epochs, the hot gas evaporates via steady winds 
from haloes smaller than a similar scale dShaviv fc Dekell 
[200^. While this is important in dwarf galaxies, it cannot 
be very relevant to the bi-modality at Mshock- 



5.2 Above the shock- heating scale 

(d) AGN feedback. The feat that AGNs exist preferen- 
tially in haloes above i\'fcrit may be due to a lower limit for 
haloe s hosting black holes iKoushiappas. Bullock fc DekeJ 
l2004l) . to starvation of AGNs by SN feedback in haloes be- 
low Merit, or to another reason. The power emitted from 
AGNs, e.g., in their radio jets, seems to be more than nec- 
essary for keeping the gas hot. Given a black-hole mass 

AfBH ~ lO^M0Vit)O 

in a galaxy of velocity dispersion V, and 
assuming that a fraction e of this mass is radiated away, the 
ratio of energies is £AGN/-Bgas ~ 7 x 10^ '^./i^o.os ^100 • For 
e > 10^^ there seems to be enough AGN energy for affecting 
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most of the halo gas. If this energy is released during rel- 
atively quiet, long phases of self-regulated AGN activity, it 
can keep the gas hot. However, black-hole physics does not 
seem to explain the charactersitic scale of bimodality. Fur- 
thermore, had the energy ratio been a measure of feedback 
effectiveness, it would have implied a decline with mass, in 
conflict with the trend of M/ L. 

The shock heating of the gas into a dilute medium is 
likeley to make it vulnerable to heating and pushing by the 
central energy source, thus providing the trigger for effective 
AGN feedback. Simulations of winds in a two-phase medium 
demonstrate that the dilute phas e is pushed away while the 
dense clouds are hardly affected et al ] |2005l) . This be- 

haviour is likely to be generic, though the mechanism by 
which the energy release d near the black hole is spread in the 
halo gas is an open issue llRuszkowski. Bruggen fc Beeelmanl 
I2004I- I Scannapieco &: Ohl2004) . If so. the feedback efficiency 
may be driven by the relative fraction of hot gas rather 
than the actual AGN energy. Figure 6 of lKeres et alJ (l2004l) 
shows that the hot fraction varies roughly oc M^''^, implying 
M/L oc M^^'^ near and above Mshock, in qualitative agree- 
ment with the observed trend item k). In this scenario, 
the shutdown scale arises naturally from the shock heating. 

(e) Two-phase Medium. Given that the cooling function 
peaks near ^ lO^K, the virialized gas at ^ lO^K develops 
a two-phase medium, with cold, de nse clouds pressure con- 
fined within the hot, dilute medium ('Field'lQea lFaU fc Ree3 
[^8^. The cloud sizes and evolution are affected by ther- 
mal c onductivity and dynamical processes llVoigt &: FabianI 
|2004^ . This can help explaining the bright-end truncation of 
the blue sequence (Mailer & Bullock 2004). Some of the gas 
is locked in the orbiting clouds and the density phot of the 
hot gas is reduced, slowing the cooling and the infall. Ap- 
proximating eq. (1171 with tcooi/icomp oc p~^T^, and recalling 
that Tv oc Mv''^, the longer cooling time makes the critical 
mass for further shock heating (fcooi/icomp ~ 1) smaller by a 
similar factor, Mshock pf/^^, namely it enhances the shock 
stability. The gas may be kept hot over longer periods by 
repeating shocks due to continuous accretion into the halo, 
which may alleviate the need for AGN feedback. Still, a nec- 
essary condition for hot gas is the initial shock heating, i.e., 
being in a halo above Mshock- 

(f ) Dynamical-friction feedback. Another heating source 
is the dynamical friction acting on galaxies in a halo 
core. The energy transferred is comparable to that re- 
quired for preventing cooling flows in cluster centres 
tEl-Zant. Kim & Kamionkowski 2004). The gas response to 
dynamical friction, unlike the DM response, has a sharp 
peak near a Mach number of unity jQstrikeilll999l. Fig. 3), 
namely when the gas is heated to near the virial tempera- 
ture in M > Mshock haloes and not in smaller haloes hosting 
cooler gas. As groups occur above a critical halo mass that 
roughly coincides with Mshock at z 1, this feedback ap- 
pears almost simultaneously with the hot medium, which 
then serves as the vulnerable victim of the same feedback 
(as for AGN feedback) . 

Figure |S| is an illustration of the strength of the differ- 
ent feedback processes, referring to the gas fraction that 
could have been heated at z ~ 0. The figure highlights 
the fact that different processes dominate below and above 
Mshock ~ 10^^ Mq. The transition from cold to hot infall has 
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Figure 8. The "strength" of the various feedback processes 
at z = 0, e.g. referring to the fraction of the initial gas that 
has been heated or removed (schematic). Different feedback pro- 
cesses are effective below and above the critical shock-heating 
scale ^ 10^^ Mq, and the feedback efficiency is at a minimum 
near this scale, giving rise to a minimum in M/ L there. 

a crucial role in determining the feedback efficiencies near 
the critical mass; it induces a minimum in feedback efficiency 
at a critical scale Mfdbk ~ Mshock and drives the shapes of 
the curves about this minimum. At higher redshifts this min- 
imum becomes wider and deeper but centered on a similar 
critical mass. 



6 THE ORIGIN OF BI-MODALITY 

6.1 A scenario from the assumed ingredients 

We propose that the cold ffows and shock heating play a 
key role in producing the observed bi-modality features. 
These features are emphasized by the similarity between the 
scales associated with shock heating, feedback and cluster- 
ing. Based on our current understanding of these physical 
processes, we assume the validity of the following: 

(a) A new mode of star formation. The collisions of the 
(partly clumpy) cold streams with each other and with the 
inner disc are assumed to produce starbursts, analogous to 
collisions of cold gaseous discs or clouds. These collisions 
are expected to produce isothermal shocks, behind which 
the rapid cooling generates dense, cold slabs where the Jeans 
mass is small. While the details are yet to be worked out, we 
assume that such a mode of star formation may be responsi- 
ble for much of the stars in the universe. In some cases, this 
can be an enhanced quiescent mode, leaving the disc intact 
without producing a big spheroid, an in other cases it may 
resemble the starbursts associated with mergers. 

(b) Hot forever. Once the gas in a massive halo is shock 
heated to near the virial temperature, it is assumed to be 
hot forever. This is based on the slow cooling time of the 
dilute hot medium and its vulnerability to AGN feedback, 
while cold, dense clouds and streams could be better shielded 
against winds and ionizing radiation. The shock heating is 
thus assumed to trigger a shutdown of all modes of star 
formation in haloes where cold streams do not prevail. 

(c) Cold streams in a hot medium. Cold streams in 
haloes above A/shock (© are assumed to supply cold gas for 
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Figure 9. Schematic illustration of the origin of the bi-modality 
in color (or SFR or stellar age) as a function of halo mass. At z > 2 
(top), continuous gas supply, regulated by supernova feedback, 
makes the galaxies evolve along the blue sequence, which extends 
beyond the shock-heating mass due to cold streams in hot media. 
At 2 < 1 (middle) , in the absence of cold streams above Mghock i 
the shock-heated gas is kept hot by AGN feedback, gas supply and 
star formation shut down, and the stellar population passively 
turns red & dead. Gas-poor mergers stretch the red sequence 
toward larger masses. When the halo mass is replaced by stellar 
mass (bottom), the red sequence is stretched toward small stellar 
masses due to satellite galaxies sharing a common halo. The color 
is correlated with the environment density via the halo mass, with 
the minimum group mass being comparable to Afshock- 



further disc growth and star formation, preferentially before 
Zcrit ~ 2 and in isolated galaxies. After Zcrit, especially in 
groups, cold streams are suppressed and a complete shut- 
down of star formation is assumed to follow. 

The bi-modality in color (or stellar age, or SFR) versus 
mass, the correlation with the environment and the evolu- 
tion with redshift, all emerge naturally from the early effi- 
cient star formation followed after 2 ~ 2 by the abrupt shut- 
down in haloes above Mghock, which typiclaly host groups. 
This is illustrated in Fig. |U] 

(a) The blue sequence is dominated by galaxies in haloes 
below Mahock, as they grow by accretion/mergers. Cold flows 
lead to early disc growth and star formation, which is regu- 
lated by SN feedback over cosmological times. Galaxies can 
get very blue because of repeating starbursts due to the 



clumpy gas supply and the interplay between infall, star- 
burst and outflow. 

(b) Bright blue extension. Some galaxies continue to be 
fed by cold streams even when they are more massive than 
Afshock, extending the blue sequence beyond L,. This occurs 
especially at z > Zait ~ 2, when the streams feeding high-cr 
haloes are relatively narrow and dense, resulting in massive 
starbursts. 

(c) The red sequence. Once a halo is more massive than 
Mshock, halo gas is shock heated; it becomes dilute and vul- 
nerable to AGN feedback. At z < Zcrit, where cold streams 
are suppressed, gas supply from the host halo shuts off, pre- 
venting any further growth of discs and star formation. If 
residual cold gas has been consumed in earlier mergers, there 
is a total shutdown of all modes of star formation above 
Afshock, allowing the stellar population to passively turns 
"red and dead" into the red sequence. The massive tip of the 
blue sequence at z > z^it becomes the massive tip of the red 
sequence at z < Zait, explaining the bright-end truncation 
of the blue sequence at low z, and the apearance of brighter, 
very red galaxies already at z ~ 1. Subsequent growth along 
the red sequence is induced by gas-poor mergers uncontam- 
inated by new blue stars. 

(d) Faint red extension. When color is plotted against 
stellar mass, the bi-modality extends to smaller galaxies 
which are typically satellites of the central galaxies in com- 
mon haloes. In haloes below Mshock, accretion onto satellite 
galaxies can keep them on the blue sequence for a while. 
In haloes above Mshock, where gas supply stops and the en- 
vironment density is high, the satellites too become red & 
dead. 

■ an associated paper dCattaneo et alJl200!TL COS) we 
have implemented the proposed new physics in a h ybrid 
SAM/N-body simulation (GallCS iHatton et aljl20o3) . The 
main new ingredient is a shutdown in gas cooling and star 
formation above Mshock ~ IO^^A/q after Zait ~ 3, while al- 
lowing for efficient star formation by cold streams in more 
massive halos prior to Zcrit- The revision yields excellent fits 
to the observed features at low and high redshifts (COS). Fig- 
ure^] shows one example of the results — color-magnitude 
diagrams which demonstrate the success of this model along 
the lines envisioned in Fig. El The top panel highlights the 
main deficiencies of the "standard" model at z = 0: an excess 
of bright blue galaxies accompanied by a sh ortage of red- 
enou gh galaxies compared to the SDSS data (|Baldrv et alJ 
|2004). The "new" model puts the blue and red sequences 
where they should be at z = 0, with a proper truncation 
at the bright-blue end and appropriately red colors in the 
red-sequence. The z — 3 diagram shows the predicted bright 
blue galaxies (which were absent in the "standard" model, 
COS). The colors distinguish between galaxies in haloes be- 
low (blue) and above 2 x lO^^M©, comparable to the mass 
separating haloes hosting field galaxies and groups. This in- 
dicates that the model recovers the correlation between color 
and environment density. It also recovers the bi-modality in 
bulge-to-disc ratio (COS). 

6.2 Origin of Observed Features 

The color-magnitude bi-modality emerges from the color- 
mass bi-modality lying at the basis of the model. The 
red satellites extend the luminosity range over which the 
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Figure 10. Color-magnitu de diagrams from a liybrid semi- 
analytic/N-body simulation iCattaneo et aIJl2005^ demonstrat- 
ing the success of the proposed scenario in reproducing the main 
bi-modality features along the lines of Fig. |5] Top: Based on 
the "standard" version of the GallCS SAM. Middle: The result 
of incorporating a shutdown in star formation above Mghock ~ 
IO^'^Mq after 2;crit ~ 3, and allowing for star formation even in 
> Msijock halos prior to Zcrit- Bottom: Same at z = 3. Blue dots 
refer to galaxies in haloes below Mghock i while red and green dots 
are galaxies in halos above Mghocki central and satellite galax- 
ies respectively. The halo mass correlates with number density in 
the environment. Contours mark number density of galaxies in 
the diagram. The color curves mark the main bodies of the blue 
and red sequences in the SDSS data iBaldrv et al.ll2'o04^ . 



blue and red sequences co-exist and highlight the gap in 
color. This gap is amplified because the galaxies making the 
transition into the red sequence once their haloes become 
> Afshock tend to be the merger remnants with big bulges 
and AGNs from the red tip of the blue sequence. 

The correlated bi-modality in bulge-to-disc ratio can be 
understood in the scenario where most of the big spheroids 
in the red sequence are passively aged galaxies that have 
grown massive stellar components already in the blue se- 
quence (COS). The transition to the red sequence is likely 
to be made by galaxies with big spheroids because (a) they 
have consumed their gas in the same mergers that produced 
the spheroids, (b) these mergers tend to occur in big haloes 
hosting groups where shock heating stops the gas supply, 
and (c) these spheroids contain massive black holes that can 
keep the gas hot. 



The strong ant i- correlation of star-formation rate (and 
blue color) with the number density in the environment 
is a natural outcome of the mass dependence. The key is 
that the minimum halo mass for groups at ~ W^^Mq is 
comparable to Mshock at jz ^ 1. The strong dependence 
of cold gas supply on host-halo mass can therefore be re- 
sponsible for the distinction between the SFR in field and 
clustered galaxies (which is thus predicted to be limited to 
late times). The galaxies dominating low-HOD haloes be- 
low Mahock enjoy cold gas supply leading to discs forming 
stars. The galaxies populating groups of subhaloes, typi- 
cally in haloes above Mshock at z ^ 1, suffered starvation 
of external gas supply and lost their internal gas in mergers, 
thus stopped forming stars and passively evolved to the red 
sequence. The faint end of the red sequence, which is pref- 
erentially present in high environment densities, is due to 
the starvation of satellite galaxies in the high-HOD haloes 
above Mahock. The model naturally predicts the seondary 
bi-modlaity scon along the red sequence (along the lines of 
Igcrhnd ct al. 2004). 

The classic morphology-environment relation, tradition- 
ally attributed to the correlation of merger rate with the 
environment, may also be viewed as a result of the cold- 
flow phenomenon. New discs are predicted to form only in 
haloes below Mahock, namely in "field" galaxies, and not in 
the group haloes above Mshock. On the other hand, the fre- 
quent major mergers in groups help build the big spheroids 
preferentially there. The abrupt shutdown of star forma- 
tion above A/shock makcs the color-environment correlation 
stronger than the morphology-environment correlation. 

The environment dependence highlights an interesting 
cross-talk between the clustering and gas processes. The dis- 
tinction between haloes hosting a single dominant galaxy 
and haloes hosting groups has traditionall y been attributed 
to ga s cooling on a dynamical time scale llRees fc OstrikeJ 
Il977l) . Our shock-stability analysis helps quantifying this 
idea. However, it seems from N-body simulations that the 
gravitational DM clu stering process has a parallel role 
jKravtsov et al.ll2004) . The HOD of subhaloes develops a 
transition from single to multiple occupancy near a compa- 
rable halo mass, associated with the current M». The HOD 
of subhaloes is similar to the HOD of galaxies deduced from 
the observed correlation function for galaxies (!^ item i). 
One can conclude that in haloes above ~ 10^^ M0 the po- 
tential wells associated with the DM subhaloes provide the 
sites for the fragmented gas collapse on the scales preferred 
by cooling, thus emphasizing this scale as the minimum scale 
for groups. The coincidence between these scales is behind 
the environment dependence of the bi-modality features. We 
have discussed above the other possible role of M* in the ap- 
pearance of cold filaments in hot haloes at high z (i^J and 
in some of the feedback processes (^SJ- 

The minimum in mean halo mass-to-light ratio M/L 
near ~ can be attributed to the maximum in 

gas supply near Mshock dictated by the shock heating above 
this mass and the associated minimum in feedback there 
(Fig. |HJ. While SN feedback gets stronger toward smaller 
haloes, AGN feedback pumps energy more effectively into 
the shock-heated medium in more massive haloes. The "fun- 
damental line" due to SN feedback b elow Mshock corre- 
sponds to M/L DC M"^/^ iDekel fc Woo. 200a'l . while the 
transition from cold to hot infall indicates M/L oc M^^^ 
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above Mahock, in the ballpark of the findings from 2dF 
item k). This settles the discrepancy between the halo mass 
function and the galaxy luminosity function both below and 
above the bi-modality scale. Our model predicts a similar 
behavior of M/L{M) at z ~ 1, to be r evealed by spec- 
troscopic surveys (such as the DEEP2/3 iMadewick et alJ 
^003, ; CoiLet al. 2004.) , which have already confirmed the 
predicted invariance of the bi-modality mass in this redshift 
range (Fig.0. 

The massive starbursts at high z, primarily due to cold 
streams, are helped by an increase in Mshock itself. At high 
z, the upper bound for SN feedback becomes ^ A/shock 
(Fig. |5J, allowing for a more efficient gas accretion at 
Mshock- If /b is 0.13 instead of 0.05 at z ~ 3, the value of 
Mshock doubles (eq. I34ll 'l. Another factor proportional to /b 
enters when translating from halo to stellar mass, yielding 
a total increase of ~ 5 in the critical stellar mass at high z. 
We thus expect a strong star-formation activity at high z in 
galaxies with stellar masses exceeding ~ 10^^ Mq. 

The global star formation history could be derived from 
the predicted SFR as a function of mass and redshift, con- 
volved with the time evolution of the halo mass function 
in the given cosmology. With the prediction that haloes 
^ Mshock are the most efficient star formers, and with the 
Press-Schechter estimate that haloes of such a mass typically 
form at z ~ 1 (Fig.|5J, the star formation density is predicted 
to peak near z ~ 1, with a relatively flat behavior toward 
higher redshifts and a sharp drop toward lower redshifts, as 
observed. In particular, the cumulative stell ar density seems 
to st op growing quite abruptly at 2 ~ 1 jPickinson et alJ 
|20o3), when the typical forming haloes become larger than 
Mshock and the cold flows are suppressed. The growth of the 
hot fraction as a function of halo mass jKeres et al J 120041) 
can be translated to the drop in star-formation rate after 
z ^ 1. The "downsizing" of star formation in galaxies is 
also helped by the shutdown of star formation above A^shock 
while smaller galaxies can make stars also later, until they 
fall into bigger haloes. 



7 DISCUSSION: OTHER IMPLICATIONS 

Possible implications on open issues in galaxy formation 
where further study is desired are worth mentioning. 
X-rays. The shock-heating scale may be detectable in soft 
X-rays, as the lower limit for galaxies and groups contain- 
ing hot halo gas. In turn, the suppression of heating be- 
low Mshock explains the m issing soft X-ray background jPenl 
ll999l:lBenson et al.ll200(]l) . We predict a noticeable suppres- 
sion of X-ray emission in the range 5 x 10^ to 2 x lO^K. 
Lyman-alpha emission. The cold (~ 10'*K) flows may 
instead be an efficient source of Ly-a radiation, possi- 
bly a ssociated wi th obse rved Ly-g emitters at high red- 
shift llKurk et alj [2003: N ilsson et all IpOOSfl . It ha s been 
argued based on SPH simulations iFardal et al.l l200lt 
IPurlanc tto ct al. 2003) that the flows radiate their infall en- 
ergy mostly in Ly-a before they blend smoothly into the 
discs. Our Eulerian simulations (ZBDK05) indicate that the 
streams do eventually shock in the inner halo. This produces 
X-ray, but, given the high density there, the X-ray radiation 
is likely to be confined to an ionized Stromgren sphere of a 
few kiloparsecs. This energy eventually transforms into Ly- 



a radiation, which could propagate out via thermal broad- 
ening and systematic redshifts. A study involving radiative 
transfer is required. 

Damped Lyman-alpha systems. The possible associ- 
ation of massive cold fiows with damped Ly-a systems 
JProchaska et al.ll20oil should be addressed in cosmologi- 
cal simulations. 

LIRGS. Cold fiows may expla in the massive starbu rsts as- 
sociated with LIRGs at 2 ^ 1 iHammer et al.ll20o3) . If half 
the stars in today's discs were formed in such LIRGs, and 
the majority of galaxies today are fragile discs, than many 
of the LIRGs could not have been produced by violent ma- 
jor mergers. The cold streams may provide a less violent 
starburst mechanism not associated with the destruction of 
discs. Simulations that properly incorporate cold streams 
should be confronted with these data. 

Angular momentum. The proposed scenario may set the 
stage for solving the angular-momentum puzzle — the over- 
production of low angular momentum spheroids in current 
cosmological simulations (Navarr o fc St cinmotz 2000). The 
solution should involve the removal of baryons with low 
angular momentum. In small galaxies, SN feedback can 
blow the gas away from their small building blocks, which 
are otherwise the main s ource of low angular momentum 
via many mmor mer gers jMaller. Dekel fc Somerville 20ol: 
iMaller fc DekeJbOOal) . In galaxies near Mshock, we find from 
cosmological simulations (ZBDK05, see Fig. |5J| that the 
low angular momentum gas is typically associated with the 
shock-heated medium, which can be prevented from cool- 
ing by AGN feedback. The cold streams come from larger 
distances with ~ 50% higher specific angular momentum, 
appropriate for producing extended discs. The feedback ef- 
fects, both below and above the critical scale, have not been 
properly simulated yet because of incomplete treatment of 
the micro-physics. 

Cold clouds. The formation of discs by a clumpy cold 
gas phase may have the following implications, (a) It may 
explain the starbursts responsible for very-blue galaxies, 
(b) It may help explaining the bright-end truncation of 
the luminosity function (JJSJ. (c) The dynamical friction 
bringing the clouds into the disc transfers energy into 
the halo, which may help explaining the discrepancy 
between the steep cusps predicted by N-body simulations 
and the fiat cores indicated by rotation curves in low- 
surface-brigh tness galaxies jPekel. Devor fc Hetzron il 

.200i iDekel. Arad. Devor fc BirnboimI ~ | 2003| : 

E^ant^Ioffman, Primack. Combes fc Shlosmanl l2004l : 
Ma fc Bovlan-Kolchin, .2004 ) . (d) The same process may 
lower the predicted maximum rotation velocity in discs 
at a given luminosity, balancing the adiabatic contraction 
of the dark halo, and thus repair the zero-point off set in 
models of the Tully-Fisher relation (e .g. iKlvpin et alll200^ : 
lAbadi et all2003l:lDutton et al.l20'M) . (e) This may explain 
the lack of anti-correlation between the residuals in velocit y 
and radius at a given luminosity JCourteau fc Rixlll999D . 
indicating comparable contributions of disk and dark halo 
to the gravitationa l potential at the effective disc radius 
jDutton et alJl2004^ . 

Dust lane. Edge-on discs above Vv — 120 km s^^ 
show a well-defined dust lane, while less massive 
discs show diffuse dust above an d below the disc 
jPalcanton. Yoachim fc Bernsteirl2004^ . A thick, turbulent. 
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Figure 11. Shock stability as a function of mass at 
Bottom: The ratio of rates, icool/*comp, as a function of halo 
mass derived at the disc (Z = 0.1, Ms = 0) and at the virial radius 
(Z = 0.03, Us = 0.15). The cooling rate is assumed to vanish for 
T < lO'^K. A stable shock is possible once icool/icomp > 1. Top: 
An illustration of the evolution of shock radius between the disc 
vicinity and the virial radius as the halo mass grows in time, based 
on the stability criterion shown in the bottom panel. 



dusty gas phase is indeed expected when SN feedback is ef- 
fective, and when cold streams shock and produce stars, 
both predicted below an appropriate scale. 
Shock heating in dwarf haloes. The cold/hot infall 
and feedback processes are expected to give rise to two 
scales characterizing dwarf galaxies. The lower bound at 
K ~ 10 - ISkms"^ (e.g., iDekel fc Woa ,2003. Fig. 3) is 
commonly attributed to the drop in atomic cooling rate be- 
low ~ 10'»K. We propose that this involves shock heating, 
in analogy to Mshock discussed above. Figure 1111 (bottom) 
shows the quantity relevant to shock stability, tcooi/icomp 
[eq. 1171 1. as a function of halo mass, now stretched to low 
masses. The molecular-hydrogen cooling rate, relevant be- 
low 10*K, is weaker and may actually be eliminated after 
g ^ 10 due to molec ule dissociation by the UV background 
iHaiman et alj[l996ll . The stability is evaluated at 2; = 
both in the disc vicinity (r — O.lRv, Z = 0.1, Us ~ 0, and 
near the virial radius {Z — 0.03, -its = 1/7). 

The top panel illustrates the deduced evolution of shock 
radius between the disc and the virial radius as the halo 
mass grows in time. Haloes below ~ 10^ M© have stable 
virial shocks. As the halo grows to ~ 6 x 10* M©, the con- 
ditions near the disc become unfavorable for a stable shock, 
but the virial shock persists. Only when the halo becomes 
~ 6 X IO^Mq shock heating stops completely. Thus, the 
shock heating prevents star formation in smaller haloes even 
before the epoch of re- ionization. Once the first stars form in 
the dense centres of big clouds, their UV radiation reionizes 
the gas, dissociates molecules, and helps preventing cool- 
ing. This provides a natural explanation for the absence 
of luminous galaxies with haloes below ~ 10^ Mq, predict- 



ing a large population of mini haloes with Vv < lOkms"^ 
which are completely dark. There is only a narrow window, 
6 X 10" — 6 X 10^^ Mq , for haloes that allow cold fiows and can 
form luminous discs at low redshifts. This range is somewhat 
broader at high redshifts, with the lower bound dropping 
below 10^ Mq and the upper bound rising above 10^^ Mq 
prior to 2: ~ 3. Within this window, it seems that below 
Vv ~ 30kms^^ most of the haloes are totally dark® and the 
others are populated by gas-poor dwarf spheroidals. This is 
the scale predicted by photoionization feedback (JU- 



8 CONCLUSION 

8.1 Summary of results 

The classic argument of cooling on a dynam ical time scale 
jRees fc OstrikeJll977tfWhite fc Reejll978D . with order-of- 
magnitude estimates of the time scales involved, provided 
an inspiring qualitative upper bound for luminous galaxies, 
at a halo mass of M ~ IO^^'^^ Mq. An analytic study o f the 
actual shock-heating process (iBirnboim fc Deke]||200l and 
this paper) now yields a more concrete halo critical scale 
at M ~ 6 X 10^^ Mq, somewhat smaller than the original 
estimate. The criterion for critical shock stability. 



cool •'comp 



(42) 



is a balance between the cooling rate and the post-shock 
compression rate, which restores the pressure supporting 
the shock against gravitational collapse. The compression 
time is somewhat larger than the crossing time at the shock 
position. The absolute magnitudes of these time scales are 
irrelevant - they could in principle both be longer than the 
Hubble time, because what matters for shock heating versus 
cold flows is only the relative rates of the competing pro- 
cesses. The most relevant critical scale is obtained in the in- 
ner halo, because as the halo grows, the shock first becomes 
stable in the inner halo, and it then propagates outward to 
the virial radius. Haloes of mass below the threshold mass 
build discs in their centres by cold flows, while in haloes 
above the threshold much of the gas is shock-heated. These 
results are confirmed by spherical hydrodynamical simula- 
tions. The same phenomenon is seen at a comparable scale 
in general cosmological hydrodynamical simulations. They 
reveal that in haloes near the critical scale, and even in 
larger haloes preferentially at z ^ 2 and especially in field 
galaxies, cold streams along the filaments feeding the galaxy 
penetrate through the hot medium, and allow further disc 
growth and star formation. 

The interplay between these cold fiows and shock- 
heating, the gravitational clustering scale, and the differ- 
ent feedback processes acting below and above a similar 
mass scale, is proposed to be responsible for the robust bi- 
modality imprinted on the observed galaxy properties. Cold 
fiows in haloes much bigger than the clustering scale allow 



® This can be deduced from the discrepancy between the flat 
faint-end luminosity function and the steep halo mass func- 
tion predicted in ACDM, given the TuUy- Fisher like velocity- 
luminosity relation of dwarf galaxies (e.g., iDekel &: Wodl2003l 
Fig. 3). 
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massive starbursts at z ^ 2, while shock heating in com- 
parable haloes at later times shuts off star formation and 
leads to big red galaxies. While supernova and radiative 
feedbacks regulate star formation below the critical scale, 
the presence of dilute, shock-heated gas in more massive 
haloes allows the AGN feedback (or another energetic source 
capable of affecting big galaxies) to keep the shock-heated 
gas hot and prevent further disc growth and star formation. 
The observed bi-modality and many of the related phenom- 
ena are argued to arise naturally from such a scenario ([0. 
The shock-heating process also plays a role in introducing 
a lower bound to haloes hosting galaxies, at ~ Mq. The 
mass range where disc galaxies can form today turns out 
to be quite narrow, between a few times 10^ to slightly 
below IO^^Mq. 

8.2 Re-engineering of SAMs 

Once the new physical processes are properly incorporated 
in the detailed models of galaxy formation, they solve many 
of the apparent conflicts between theory and observation. 
At a crude level, one might have naively thought that since 
the cooling time is anyway assumed to be shorter in smaller 
haloes, the details of the cold-ffows and shock heating would 
not matter much to the final result. However, a closer inspec- 
tion reveals that there are several key features which make 
a qualitative difference: 

(a) Star formation. The near-supersonic cold streams pro- 
vide a new efficient mechanism for early star formation. This 
is in contrast to the gradual infall of cooling shock-heated gas 
assumed in most SAMs, which starts from near rest, accretes 
smoothly into the disc, and joins the quiescent mode of star 
formation there. We find that the expected cold-gas supply 
is significantly more efficient than assumed in most current 
models even in small haloes (Cattaneo, Neistein, Birnboim 
& Dekel, in preparation). 

(b) Heating inside out. The concept of an expanding 
"cooling radius" used in current SAMS is limited to mas- 
sive haloes where a virial shock exists. Otherwise, it is the 
shock causing the heating which is propagating from the 
inside out. 

(c) Shutdown of star formation. The combination of 
shock heating and AGN feedback provides a mechanism for 
shutting off disc growth and star formation above a thresh- 
old halo mass. 

(d) Cold streams. Cold streams that penetrate through the 
hot media continue to make discs and produce stars in haloes 
above Mghock- This happens mostly at z ^ 2, and prefer- 
entially in less grouped galaxies, allowing big blue galaxies 
mostly at high z and some at low-density environments, and 
enforcing a sharp shutdown of star formation at late times 
and especially in clustered galaxies. 

A practical schematic recipe for the critical halo mass 
below which cold streams prevail and above which one may 
assume a shutdown of gas supply and star formation is: 

{A-^shock, Z < Zcrit 
A/, Jm^\ (43) 

where the critical redshift Zcrit is defined by /M*(2crit) = 
Mshock, the clustering scale Mt{z) is given by eq. (IA18II . and 
/ is a numerical factor of order a few. Our best estimates 



for the parameters are log Ms ~ 11.8 (but possibly another 
value in the range 11.3-12.3) and / — 3 (possible range: 1- 
10). Using the approximation logM, ~ 13.1 — 1.32: (2 2), 
we obtain 2crit — 1.4 for / = 3. This recipe should allow big 
blue systems at 2 > 2, eliminate big blue systems and make 
big red galaxies at 2 ^ 1, and generate a bi-modality near 
Afshock. This scheme can be refined to allow for a smooth 
transition above the critical scale by applying the shut-down 
to a varying fraction of the gas and by breaking the streams 
into clumps which will generate high peaks of starbursts. 

In addition, one may wish to have an effective minimum 
requirement for the central black-hole mass in order to en- 
sure enough feedback energy for maintaining the gas hot. 
This may emphasize the bi-modality gap in color and bulge- 
to-disc ratio. However, the proposed shut-down by halo mass 
may be enough for ensuring sufficient bulge mass and black- 
hole mass. 

The SAMs should be re-engineered to incorporate these 
processes and thus help working out the detailed impli- 
cations of the proposed scenario. Preliminary attempts to 
do that, using two different SAMs, indicate that the in- 
corporation of the new proposed processes outlined above 
indeed leads to significantly better fits with the observed 
bi-modality features along the lines proposed in SJHKCOS). 

8.3 Open issues 

In parallel, the physics of the involved ingredients should 
be properly worked out in more detail, starting with the 
following two hypotheses that were laid out in ^ 

(a) Fate of cold streams. A detailed investigation is re- 
quired of the way by which the cold streams evolve and even- 
tually merge with the central disc, the associated star for- 
mation, and the resulting feedback process. While progress 
can be made using toy models and simplified simulations, a 
proper analysis will require simulations of higher resolution 
than are currently available. In particular, whether or not 
the predicted star bursts could be associated with the big 
dusty sources indicating massive sta r formation at high red - 
shifts, such as the SCUBA sources llChapman et alJl2003H . 
remains to be determined once the theory is worked out and 
the observed characteristics of these sources are clarified. 

(b) AGN feedback. The physics of AGN feedback is an- 
other unknown. One wishes to understand how the available 
energy originating near the central black hole is transferred 
to the hot gas spread over the halo. The physics of how ther- 
mal conductivity may heat the gas is also to be investigated. 
The increased efficiency of these feedback mechanisms in the 
presence of a hot medium as opposed to their effect on cold 
flows and clumps are to be quantified. 

Parallel attempts to work out the details of the physical 
input and to incorporate it in quantitative models of galaxy 
formation will lead to progress in our understanding of the 
galaxy bi-modality and the associated features. 
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APPENDIX A: USEFUL RELATIONS 

We summarize here the cosmological relations used 
in the analysis of ^ T his is rather basic materia l, 
based for example o n Lahav. Rees. Lilie fc Primack' ("1991"); 
ICarroU. Press fc Turner. (il992h and lMo fc White (,,2002) . By 
specifying it here in a concise and convenient form, we hope 
to allow the reader to reproduce our results and use them in 
future work. Additional relations associated with the spher- 
ical top-hat collapse model are brought in the appendix of 
BD03. 



Al Cosmology 

The basic parameters characterizing a flat cosmological 
model in the matter era are the current values of the mean 
mass density parameter fim and the Hubble constant Hq. 
At the time associated with expansion factor a = 1/(1 + 2:), 
the vacuum-energy density parameter is flA{a.) = 1 — f2m(a) 
and 

= o "ro"'-3 ■ (Al) 

The Hubble constant is 

H{a) = Ho {Qa + Qma-'^)^^", (A2) 
and the age of the universe is 

tu„iv(a) = -i/(a) (|i_n^(„)|)i/2 • (A3) 

The mean mass density is 

Pu ^ l.SSxlO'^^ Qh^a'^ ~ 2.76xlO'^"ano.3 feo.Td"^ , (A4) 

where Slmo.3 = fim/0.3, h = Ho/100 km s^^ Mpc~^ , and 
ho.T = h/0.7. 

A2 Virial relations 

The virial relations between halo mass, velocity and radius, 



2 GMv Mv 



(A5) 



become 
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Mil ==; 6.06 ViloA^^^ - 342 (A6) 

where Mn = A/v/IO^Mq, Vioo = K/100kms~\ Ri = 
7?v/l Mpc, and 

A = (A200 fimo.a /io.t)"'^^ a- (A7) 

An approximation for A (a) in a flat universe (Bryan & Nor- 
man 1998) is: 

A(a) ~ (IStt^ - 82!:!A(a) - 39QA{af)/^lmia) . (AS) 
The virial temperature can be defined by 

For an isotropic, isothermal sphere, this equals a^, where a 
is the one-dimensional velocity dispersion and the internal 
energy per unit mass is e = (3/2)cr^. Thus 

Vfoo-2.79r6 Mil ~ 28.2 Tg/^A^/^ (AlO) 
where Ts = n/lO^K. 



where M and the comoving radius R are related via the 
universal density today: M = ^poR'^- The mass of 2-a 
fluctuations is obtained by se tting v(M, a) = 2, etc. B ased 
on the improved formalism of ISheth fc TormenI ||20o3), the 
fraction of total mass in haloes of masses exceeding M is 




(A19) 



This fraction for 1-a, 2-a, and S-cr fluctuations is 22%, 4.7%, 
and 0.54% respectively. 

Figure|5|shows the PS mass M* as a function of redshift. 
For the standard ACDM with erg = 0.9 its value at z = 
is M*o = 1.36 X 1O^^M0. One can see that an excellent 
practical fit in the range ^ 2 ^ 2 is provided by a power 
law in this semi-log plot: logM* « 13.134 — 1.32. At larger 
redshifts this gradually becomes an underestimate. Trying 
to provide crude power-law approximations, we find that 
M, oc a*'^ cx t^'^ are crude approximations in the range 
^ 2 ^ 1, and that Af, cx a"' cx t'^ are good to within a 
factor of 2 in the range ^ 2 ^ 2. These power laws become 
overestimates at higher redshifts. 



A3 Press Schechter 



Linear fluctuation growth is given by llLahav et al.lll99lt 
ICarroU et al]ll992t IMo fc Whit^l2002h 



N q(a,) 



where 
9ia) 



^f^m(a) 

o ^o^-i/T „ (l + nm(a)/2) 



(All) 



(A12) 



(l + J7A(a)/70) 

The CDM po wer spectrum is approximated by 
jBardeen et alJll986h : 



P{k) cx kT^{k), 
with 

T(fc) = 1^(1 + 2.34,) 



(A13) 
(A14) 



2.34 g 

X [1 + 3.89g + il6.1qf + (5.46g)^ + (6.71g)'']"'''*, 

where 

q = k/{Q.^h^Mpc'^). (A15) 
It is normalized by as at i? = 8 ^~^Mpc, where 



'^^{R) = ^l dkk^ P(k)W'\kR), 
27r 



(A16) 



and with the Fourier transform of the top-hat window func- 
tion 



iy(a;) = 3(sin x ~ x cos x)/x 



(A17) 



In the Press Schechter (PS) approximation, the charac- 
teristic halo mass M,{a) is defined as the mass of the 1-a 
fluctuation, 



1 = u{M, a) = 



Sc 



D{a)a{M)' " 
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1.69, 



(A18) 



